Monopole core instability and Alice rings in spinor Bose-Einstein condensates 
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We show how the length scale hierarchy, resulting from different interaction strengths in an 
optically-trapped spin-1 ^^Na Bose-Einstein condensate, can lead to intriguing core deformations 
in singular topological defects. In particular, a point defect can be unstable with respect to the 
formation of a stable half-quantum vortex ring (an 'Alice ring'), providing a realistic scheme to use 
dissipation as a sophisticated state engineering tool. We compute the threshold for stability of the 
point monopole, which is beyond the current experimental regime. 
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The rich order parameter space of multi-component 
Bose-Einstein condensates (BECs) can admit truly 3D 
topological excitations 0, S S S 0^ beyond the sim- 
ple quantized vortices of single-component BECs. Such 
structures are of interest in a wide range of physical 
contexts, but dilute atomic BECs offer the unusual ad- 
vantage that we can fully explore, e.g., the short-range 
physics in topological defect cores, where the order pa- 
rameter may explore a larger space than the usual ground 
state manifold. In this Letter we show how this can result 
in rich and surprising core structures, by demonstrating 
a spontaneous deformation of a singular point defect to 
an energetically stable half-quantum vortex ring. 

We examine the recently presented case j5| of a defect 
analogous to the 't Hooft-Polyakov monopole , in a an- 
tiferromagnetic, or polar, spin-1 BEC 7]. We will show 
that it is only in the strongly antiferromagnetic regime, 
which is not attained in current experiments, that its 
stable core will be the spherically symmetric hedgehog, 
with a total density depression, of Ref. 5]. In the weakly 
antiferromagnetic regime that currently holds, the wave- 
length at which the antiferromagnetic constraint may be 
violated is much larger than that at which the total den- 
sity constraint fails. The stable defect core therefore ex- 
tends to this larger size, and holds non-zero average spin 
instead of a density zero. The singular point defect itself 
deforms to a circle: a half-quantum vortex ring (Figs. Q 
and PI , called an 'Alice ring' by high energy physicists 
H, id linj , which carries a topological charge similar to 
delocalized magnetic 'Cheshire' charge This forms 
an interesting connection between ultra-cold atom exper- 
iments and elementary particle physics. It also shows 
that dissipation, often an obstacle in state engineering, 
can sometimes perform the intricate final step in produc- 
ing an exotic object. 

We consider the BEC of spin-1 atoms. In the absence 
of a magnetic trapping potential, the macroscopic BEC 
wave function is determined by a spinor wave function ^ 
with three complex components 
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FIG. 1: The stable half-quantum vortex ring (Alice ring), 
when the energy of an initial spherically symmetric monopole 
was minimized by continuously deforming the field configu- 
ration. The asymptotic distribution of the spin quantization 
axis d(r) (left), for r ^ ^a, forms the radial hedgehog. For 
visualization purposes, the unoriented d(r) field is drawn by 
cones. We show the constant surface density plots (right) for 
|'0i(r)|^ (red) and for |'0_i(r)|^ (blue), where the monopole 
core is deformed for r < with the two line vortices sepa- 
rating. 



field theory for this system reads: 
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where F are the 3-by-3 Pauli spin matrices, (F) = 
. F . denotes the average spin, p = the to- 
tal atom density, cq = 47rfi^(2a2 + ao)/3m, and C2 = 
47rfi^(a2 — ao)/3m, for ap the two-body s-wave scatter- 
ing length in the total spin F channel 12]. For ^^Na, 
aa {a2 - ao)/3 2aB and (2a2 + ao)/3 SOa^, 
indicating C2/C0 :^ 0.04, where as = 0.0529nm is the 
Bohr radius. Here V denotes the external potential, for 
an isotropic optical dipole trap with the frequency co: 
V(r) = muP'T^ 12. For ^^Na C2 > 0, and the energy is 
minimized by setting (F) = throughout the BEC for 
the case of a uniform order parameter field. The zero av- 
erage spin corresponds to the ground polar state, where 
we may determine all the degenerate states by means of 
the macroscopic BEC phase and a real unit vector d(r) 
defining the quantization axis of the spin. The BEC wave 
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FIG. 2: The spin profile of the Alice ring displayed in Fig. Q] 
The spin expectation value (F) (left) is nonvanishing along 
the half-quantum vortex ring core. The absolute value of the 
spin |(F)|^ (right) between the isosurface sections is indicated 
by a colormap from red (|(F)|^ = 0.35) to purple (|(F)|^ = 1). 
To display the isocaps of |(F)|^ we have rotated the surface 
with respect to the spin field graph (F) approximately 7r/4 
counter-clockwise along the z axis and cut the ring half. 



function then reads: 
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As in the similar polar phase of superfluid ^He-A 0, 
[15I H^ . however, the states (d, (/?) and (— d, + tt) are 
identical The polar order parameter space, which 

may appear to be 5*^ x S*^, is actually factorized by the 
two-element discrete group Z2. Consequently, we take 
the d field to define unoriented axes rather than vectors. 

The spherically symmetric monopole is obtained 
from Eq. Q by the radial hedgehog field d(r) = r = 
(sin ^ cos 0, sin ^ sin 0, cos ^) with = and p = 

puir^ minimizing the energy of the symmetric config- 
uration 5]. This is singular at the origin indicating a 
point defect with p(0) = 0. Spinor component ?/^o resem- 
bles a dark soliton and form perfectly overlapping, 
straight singly-quantized vortex lines with opposite cir- 
culation, perpendicular to the phase kink plane. The 
topologically invariant winding number. 
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is defined on any closed surface dVt that encloses the ori- 
gin. Because the sign of d is ambiguous, though, the 
sign of W is arbitrary. Moreover, because the d field is 
actually unoriented, the monopole point defect may be 
continuously deformed into a circular line defect, without 
changing W (on any surface enclosing the singular ring), 
by punching a hole in the spherically symmetric core; 
see Fig. |2l To keep ^ single- valued on the disc bounded 
by the ring, the macroscopic phase must change by tt 
around any loop that links the defect circle, while there 





FIG. 3: Continuous deformation of the radial hedgehog into 
an Alice ring. Lines represent flow lines of the order parameter 
d field. We show planar sections of the monopole unperturbed 
(left) and deformed into an AUce ring (right). Note that the 
asymptotic behavior of d remains unchanged. 



is also a 7r-disclination in d on the disc. We identify this 
structure as a half-quantum vortex line , forming a 

closed circular ring, also called an Alice ring |§. 1^. HoL HH . 
To keep ^ single- valued on the ring itself, one can either 
have p vanish there, or have | (F) | = 1 on the ring instead. 

Thus, because of the ambiguity of the direction of d in 
the polar state, different core structures can be smoothly 
deformed into each other. The core stability is there- 
fore solely determined by energetic considerations. Phys- 
ically, the p = ring evidently has higher energy than 
the p = point; but energetic stability of the point 
defect against the Alice ring with |(F)| = 1 core is a 

nontrivial question. We may define two healing lengths 
= (8to,p)-V2 ^^^^^ = (g^^^^)-i/2^ i^Eq^ They 

describe the length scales over which p and |(F)|, re- 
spectively, tend to their bulk value when subjected to 
localized perturbations. (We can also say that ^5 and 

correspond to the excitation wavelengths below which 
density and spin fluctuations cease to be energetically 
suppressed.) Since at the defect core we may have ei- 
ther p = or 1(F) I = 1, and determine the core 
sizes in the two cases. We numerically demonstrate, by 
decreasing the ratio a^/as, that the spherically symmet- 
ric monopole core with the total density suppression be- 
comes unstable to formation of an energetically stable 
Alice ring with | (F) | = 1 core. We find the radial hedge- 
hog to be unstable in a linear stability analysis in bulk 
at C2/C0 ^ 0.17 and in the mean- field theory in a trap 
at C2/C0 ^ 0.16, using an experimentally feasible set of 
parameters for ^^Na. 

Before turning to the full nonlinear mean-field theory 
of the trapped BEC, we study the linear stability of the 
radial hedgehog in the homogeneous case. We expand 
the Hamiltonian with V = 0, to second order around 
^ M' Perturbations of the form 
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decouple from all others, for real Cj, rj and ^, if 7^ and 
satisfy (here X = C2/C0 and p = PM{r)/po, where po 
denotes the constant asymptotic value of the density) 
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for some real eigenvalue A. For small C = |C|, the 
change in free energy due to this perturbation in ^ will 
be SE = XC'^m/h^. Numerical solution shows 1^ that 
there exists one negative A, hence three degenerate in- 
stabilities proportional to when C2/C0 ^ 0.17. We 
can show that no other instabilities exist in bulk, so for 
C2/C0 ^ 0.17 the point defect is stable. Since the pertur- 
bation gives 
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where (F) is with respect to the state ^ m + (^^, we have 
1(F) I = 1 on a circle in the plane through the origin 
perpendicular to C, at radius r*, such that puir^) = 
[Crj{r^)Y . The form of pm{t) (rising monotonically from 
to po) and r] (decreasing monotonically to zero as 
r 00) ensures that grows monotonically with C. As 
a result, a randomly oriented Alice ring will form sponta- 
neously from the symmetric monopole, as C grows from 
zero by relaxation, for C2/C0 ^ 0.17. 

In the full 3D classical mean-field theory of the spin-1 
monopole in an isotropic trap, we minimized the energy 
by evolving the GP equations, 

i h— = (-^ V2 + V + Cop)* + C2P(F) • F • * , (7) 

in imaginary time. The initial state was an approximate 
spherically symmetric monopole solution with a point 
core, embedded in a Thomas- Fermi density profile. The 
integration was performed on a spatial grid of 128^ points 
using the split-step method. At every time step we nor- 
malized the wave function to fix the total atom number. 
The numerical simulations were fully 3D without impos- 
ing any symmetry on ^ as it relaxed. 

We varied the relative strength of the two interaction 
coefficients cq and C2, as well as the total atom number 
N. The dynamics depends on the dimensionless interac- 
tion strengths Cq = AiiNas/l and c'2 = AnNaa/l^ where 
/ = [h/{muj)]^^'^ . Hence, for a given atom with fixed val- 
ues of scattering lengths, the results are unchanged for 
any scaling of length and time, which does not change 
N'^uj. In optical dipole trap experiments on ^^Na BECs, 
a wide range of aspect ratios and trapping strengths have 
been demonstrated, with the typical values of I /as vary- 
ing between 500 and 4000 7]. 

For small values of C2/C0, the point core deformed 
into a half-quantum vortex ring. With the initial ra- 
dial hedgehog at the trap center we found vortex ring 



configurations as local energetic minima, with the size of 
the ring approximately determined by ^a!2^- In Figs. 
and 121 we display such results using the parameters of the 
spin-1 ^^Na with C2/C0 = 0.04 and c'q = 2 x 10"^. With 
the trapping frequency = 27r x lOHz, this corresponds 
to :^ 4 X 10^ atoms. In the initial state of the radial 
hedgehog, the components 2p±i form two oppositely cir- 
culating vortex lines with perfectly overlapping density 
profiles. Due to dissipation, the two vortices in sep- 
arated near the trap center, forming an Alice ring. Note 
that this corresponds to the growth of the perturbation 
C2 in Eq. The numerical noise in this example effec- 
tively seeded the instability C2 yielding 
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Since rr]/^pM is constant at the origin, but vanishes as 
r ^ 00, this means that near the center the two vortices 
in iIj±i separate in the xz plane. Growth of Ci involves 
a similar separation in the yz plane. One can clearly see 
the TT winding of the macroscopic phase around closed 
loops threading the Alice ring in Eq. ©. 

With some initial parameters we also observed the ra- 
dial hedgehog evolving to an Alice ring with the vortex 
cores in iIj±i expanding, but not separating, while the 
phase kink in tpo no longer maintained a total density de- 
pression at the center. This clearly represents the growth 
of the perturbation Cs in Eq. since the term propor- 
tional to ^ (r) is simply a perturbation in the macroscopic 
phase cp of Eq. ((2|). In the general case of a symmetric 
initial monopole, dissipation spontaneously breaks the 
spherical symmetry and the resulting ring has an arbi- 
trary orientation. 

It is easy to see that the Alice ring in bulk must sta- 
bilize at some finite value of the radius, because the 
amount of energetically costly non-zero (F) rises as the 
ring grows. This was also true in a trap for small enough 
rings. However, with smaller atom numbers, Cq ^ 10^ for 
^^Na, a larger ring resulted, which was destabilized by 
the inhomogeneous density profile: the two vortex lines 
in iIj±i completely detached and left the atomic cloud in 
opposite directions. 

The half-quantum vortex ring, shown in Figs. Q and [3 
was still unstable with respect to drifting out of the 
atomic cloud as a unit if initially displaced from the 
trap center. This is because of the reduced order pa- 
rameter bending energy at lower atom densities. How- 
ever, we successfully stabilized the Alice ring by creating 
a local density minimum at the trap center, by adding 
to the harmonic trapping potential the optical poten- 
tial Vl, simulating two orthogonal blue-detuned focused 
Gaussian laser beams: Vl = Fiexp[— 2(x^ -\-y'^)/w'^] + 
V2 exp [— 2(x^ -\- z'^)/w'^]. We note that the Alice ring can 
always be made long-living on experimental time scales 
even without the stabilizing field, since the drift time 
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should increase with the BEC size, while the ring forma- 
tion time should decrease. 

We also investigated the onset of the point defect 
core instability by varying C2 for a given Cq = 2 x 10^. 
We found the symmetric monopole to be unstable for 
C2/C0 ^ 0.16, in a good agreement with the linear stabil- 
ity analysis in bulk. The precise value of the threshold 
is very difficult to determine in a trap, since the conver- 
gence in the simulations is particularly slow close to the 
threshold. In the case of very long runs, numerical noise 
also eventually displaces the monopole so that it starts 
drifting out of the atomic cloud. For smaller values of 
C2, representing unstable Alice rings which break apart 
due to the inhomogeneous density, the instability of the 
point defect occurred with much larger values of C2/C0, 
when the trap length scale became comparable to the 
polar healing length / ~ emphasizing the nontrivial 
nature of the finite size effects. 

In the experimental preparation of an Alice ring by re- 
laxation, it is helpful to note that defects and textures in 
multi-component BECs may generally be viewed as com- 
binations of vortex lines, rings, and phase kinks 0. As 
we already noted, the radial hedgehog has an especially 
simple structure of vortex lines and phase kinks which 
have been experimentally created in individual atomic 
BEC components, e.g., by means of localized laser fields, 
which are rapidly swept around the trap during a Ra- 
man transition 0, [22, 113 • The spherically symmetric 
monopole could similarly be created, by a sequence of 
Raman pulses, in a straightforward generalization of the 
pulse sequences proposed in Ref. 1]. Dissipation plays 
a crucial role in the state engineering process and will 
then perform the final step of generating a stable Alice 
ring by spontaneously breaking the spherical symmetry. 
The stable radius of the Alice ring may easily extend 
over several microns, making the core quite observable, 
possibly even without ballistic expansion. 

We can summarize our results in basic terms. In 
^^Na BECs the weakness of the antiferromagnetic energy 
{as ^ aa)^ topology, and the gradient energy of the or- 
der parameter, may conspire to favor nonvanishing spin 
values, even for > 0: It is energetically more favorable 
to violate the antiferromagnetic constraint than to force 
the superfiuid density to vanish. As a result, a monopole 
core deforms to a ring and exhibits a nonvanishing spin 
expectation value and a non-zero superfiuid density. In 
other words, the strong order parameter bending energy 
close to the singular defect mixes the polar and the ferro- 
magnetic phases of the spinor BEC, rather than forcing 
the total superfiuid density to zero at the singularity. 
Due to the length scale hierarchy, the stable size of the 
defect core will then be determined by instead of the 
much smaller ^5. And the creation of this intricate Alice 
ring structure will occur spontaneously, from the much 



simpler symmetric monopole, by relaxation alone. 

While the motion and interaction of point defects was 
analyzed in 5], the dynamics of Alice rings remains a 
question for future theoretical and experimental study. 
So does the possible meta-stability of Alice rings for C2 ^ 
O.ITcq. Topological defects in spinor BECs promise a rich 
phenomenology, which current experimental techniques 
will allow us to explore. 
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